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We suggest an experimental measurement that could detect the appearance of tau neutrinos due 
to Vr oscillations of atmospheric neutrinos by measuring the energy spectra of neutrino 

induced showers, r neutrinos deposit a large fraction of their energy in showers generated by Vr CC 
interactions and the subsequent r-lepton decay. The appearance of Vr will enhance the spectrum 
of neutrino induced showers in energy ranges corresponding to the neutrino oscillation parameters. 
A shower rate lower than the 'no oscillation' prediction is an indication for Vy. — > Vs oscillations. 
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I. INTRODUCTION 

The Super-Kamiokande experiment, a densely instru- 
mented water-Cherenkov detector of very large dimen- 
sions, confirmed ||^ earlier indications Q that the abnor- 
mal ratio of atmospheric muon to electron neutrinos can 
be interpreted best as muon to tau or sterile neutrino 
oscillations. The oscillation hypothesis is supported by 
the low ratio of muon to electron neutrino events and by 
the disappearance of muon neutrino events as a function 
of the distance to their production site. Several inde- 
pendent data sets, measuring different neutrino energy 
ranges and neutrino interaction processes, are fully con- 
sistent with oscillations in maximum mixing and Ato^ 
value of 3.5x10-3 eV^ @. 

None of the data sets, that have been currently an- 
alyzed, however, contain any signatures of tau lepton 
appearance. The low deep inelastic scattering cross 
section below 10 GeV and the very short r-lepton life- 
time prevent the direct observation, although the exper- 
imental sample should contain a number of r-leptons. 
Several long baseline accelerator experiments have been 
proposed with the aim to detect r neutrinos and confirm 
the oscillation hypothesis. These employ an 

accelerator neutrino beam that is more restricted in en- 
ergy than the cosmic ray beam and high resolution near 
and distant ('^700 km) detectors that are able to identify 
Vt interactions. 

We propose an experiment that uses the short lifetime 
of the T-leptons to detect one of the signatures of their 
appearance. Charge current Vt interactions and the sub- 
sequent T decays will create hadronic/electromagnetic 
showers that will practically coincide in vertex and in 
time. Much higher fraction of the energy will be de- 
posited in the form of showers than in either v^^ CC in- 
teractions or in neutral current (NC) interactions of any 
neutrino flavor. A measurement of the energy spectrum 
of shower events initiated by atmospheric neutrinos will 
be able to register the appearance of Vt, that are present 
at a very low level in the atmospheric neutrino flux in the 
absence of neutrino oscillations. Q The shower signal is 
generated by atmospheric neutrinos of energy above 10 
GeV. The neutrino induced showers can be contained in- 



side a big water (or ice) Cherenkov detector that does 
not have to be very densely instrumented. 

In the case of — > v-^ oscillations an excess of neutrino 
induced showers will created at certain shower energies 
that will reflect the values of the oscillation parameters. 



In the case of Vu 



neutrinos the shower rate will 



correspondingly decrease, although by smaller amounts. 



II. SHOWERS GENERATED BY r NEUTRINO 
CC INTERACTIONS 

We envision a detector that consists of three or four 
strings of photomultipliers in the fashion of the high en- 
ergy neutrino telescopes j|] on a circle of radius 10 to 15 
meters instrumented with photomultipliers every 4 to 5 
m. Such an arrangement has been already proposed 0| 
for the measurement of the energy spectrum of neutrino 
induced muons. This detector will not be able to recon- 
struct the shower development or differentiate between 
purely electromagnetic and hadronic showers, but should 
be able to contain showers of energy up to 1000 GeV. The 
depth of maximum ^max for 1000 GeV electromagnetic 
showers is '^290 g/cm"^ and about 100 g/cm^ larger for 
hadronic showers, i.e. 3-4 meters of water or ice. Ac- 
counting for the absorption length of Cherenkov light, 
about 25 m, the electromagnetic and hadronic showers 
will appear indistinguishable in such relatively crude (in 
comparison to Super-Kamiokande) detector. So we de- 
fine as shower energy Esh the total energy released in the 
form of hadrons, photons and electrons in the final state, 
i.e. Esh = y X in NC interactions. 

Assuming that muon neutrinos and antineutrinos os- 
cillate into Vr (Dt) there will be three sources of neutrino 
induced showers: those due to Veii^e) CC interactions, 
NC interactions of all three neutrino flavors and CC in- 
teractions of Vr and Vr- Electron neutrino CC interac- 
tions deposit the total neutrino energy in the form of a 
shower. The fraction of shower energy in NC interactions 
is defined by the differential cross section da/dy which is 
energy dependent. To determine Vr{vr) CC interaction 
contribution one has to add to da/dy the fraction of the 
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T-lepton energy (I — y) x which is not carried away 
by neutrinos at its decay. 
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FIG. 1. Fraction of the neutrino energy that converts to 
shower energy for Ei, — 10 GeV. The solid histogram is for 
i^T CC interactions, dotted one is for NC interactions. 

Fig. shows the distribution of the fraction of neu- 
trino energy deposited in shower form for NC neutrino 
interactions and in CC Vr interactions calculated using 
the GRV structure functions |^ and r decays performed 
by JETSET74 § for = 10 GeV. The polarization of 
T leptons is not taken into account. The two distribu- 
tions are distinctly different. Egh/Ei, peaks at 95% in 
Vr CC interactions with an average of 66%. NC interac- 
tions yield on the average 48% of the neutrino energy in 
showers. 



III. NEUTRINO FLUXES FOR DIFFERENT 
OSCILLATION PARAMETERS 

In the energy range below 100 GeV atmospheric neu- 
trinos are generated predominantly by the decay chain 
tt''^ v^{Vf^) -t- /X ^ ^^{v^) -t- Ve{ve) + 6=*= of pions pro- 
duced in interactions of cosmic rays in the atmosphere. 
In the 100 GeV range the spectrum of neutrinos from 
pion decay is one power of E steeper than the cosmic ray 
spectrum because of time dilation, while those of neutri- 
nos from muon decay are steeper by two powers of the 
energy. 

The oscillation probability in a simple two neutrino 
scenario is given in convenient units as 



= sin^ 26lsin^ 



1.27 



(L/km)(ATOVeV2) 



[E^/GeY) 



between the two mass eigenstates. 



is the mixing angle 
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FIG. 2. Neutrino fluxes of all three flavors (dots: fe, 
dashes: u^, dash-dash: Vt) in the presence of Vf^ — > Vt 
cillations with maximum mixing and Am^ = lO"'^, 10 
and 10"^ eV^ from left to right. 
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Fig. H shows the fluxes of atmospheric v+D for all three 
neutrino flavors coming from the lowest tt stcradian of 
solid angle (cosf? from -1 to -0.5) for maximum mixing 
and Am^ — 10^^, 10~^-^ and 10~^ eV^ derived from 
the atmospheric neutrino fluxes of Ref . [0 . The flux of 
muon neutrinos is now split between muon and tau neu- 
trinos, the valleys in muon neutrinos corresponding to 
peaks in tau neutrinos. Because the neutrino pathlength 
in this cos range varies only between and 2Rq and 
the neutrino energy of interest is relatively high, the os- 
cillation patterns are not completely smeared, as they are 
in the GeV range. The peak in the Vt spectrum at 70 
GeV and Arn^ — 10^^ eV^, for example, corresponds to 
Vfj, oscillation probability of 0.83 for 1 R^ and 0.55 for 
2 R0. The energy spectrum of Ve is very steep and the 
Vt flux (that derives from v^) is higher than the flux 
below 500 (50) GeV for Am^ = 10"^ (j^g-S) ^y2^ ^^^^le 
angular range were narrower and closer to vertical direc- 
tion the peaks and valleys in the I'r flux would be even 
more obvious. 



IV. THE ENERGY SPECTRUM OF NEUTRINO 
INDUCED SHOWERS 

Fig. H shows the energy spectra of showers generated 
by different neutrino flavors and interactions. The solid 
line indicates showers generated by NC interactions of all 
three flavors and the dotted line gives the showers of Vg 
and De CC interactions. The sum of these two contri- 
butions is the expected shower energy spectrum in the 
absence of neutrino oscillations. The other three curves 
show the contribution of and Dr CC interactions for 
the three values of Arn^ used in Fig. |^. We assume here 
that the muons generated in i/^ CC interactions will be 
detected and used to veto the accompaning showers. 
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FIG. 3. Shower energy spectra for: NC interactions 
(solid) , Us CC events (dots) and fr CC events for oscillations 
with maximum mixing and Am^ = 10"^ (dashes), 10"^'^ 
(dash-dot) and 10"^ (dash-dash) eV^. 

At the lowest shower energy shown, 5 GeV, the contri- 
bution of electron neutrino CC interaction is the biggest 
because Ve deposit all of their energy in the form of show- 
ers. Because of the steepness of the v,, energy spectrum, 
however, NC interactions dominate for shower energies 
above 15 GeV. The contribution of i^r(i^r) CC interac- 
tions is 5-10% at 5 GeV. It increases with energy to be- 
come dominant for Am^ > 10~^'^ and vanishes at high 
shower energy where the neutrino oscillation probability 
is very small. The low contribution at threshold is re- 
lated to the small i^r CC cross section at energies only 
slightly higher than the r lepton mass. 
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FIG. 4. Shower energy spectra for v-r (thick lines) 

and v^jt — > Va (thin hnes) oscillations with maximum mixing 
compared to the no oscillation case (solid line). Am^ =0.01 
eV^ is shown with dashes, 0.005 eV^ - with dash-dot and 
0.002 eV^ - with dash-dash). 



Fig. ^ shows the total shower rate (the sum of all three 
contributions in Fig. % for Am^ of 0.002, 0.005 and 0.01 
eV^. The no oscillation case is given with a solid line. 
The upper three histograms (thick lines) are for i/^ — > v-^ 
oscillations while the lower three (thin lines) are for v^^ 
Vs case. The shapes of the energy spectrum that reflect 
the Vr contribution are as essential as the increase in the 
total rate. With the decrease of Am^ the shower excess 
moves to lower energy. The biggest excess for Am^ = 
0.01 eV^ is in the range of 40 to 150 GeV. Ara^ = 0.005 
(0.002) cV^ causes an excess above 20 (10) GeV. 

In the case of v-^ — > Vg oscillations the changes in the 
shower rate are less spectacular. The only difference with 
the no oscillationscase is the missing contribution of NC 
interactions of the oscillated muon neutrinos. Still, the 
decrease of the rates in the 10 - 50 GeV range is of order 
30 - 40 %. 

Table || gives the total rates of shower events in four 
wide logarithmically equally spaced energy intervals. The 
total rates above E^k = 10 GeV are 2.1 events per Kt.yr 
in TT steradian in the absence of oscillations. This rate 
increases by 60% to 3.35 for Am^ = 10~^ eV^ in the case 
of Vfj, — > oscillations and has intermediate values for 
lower Am^. The corresponding decrease in the case of 
sterile neutrinos is by 25% to 1.60. 

TABLE I. Shower event rates (in (Kt.yr)"^) for different 
Am^ in four energy bins for the lowest tt steradian of solid an- 
gle, i.e. between the nadir and 30 degrees below the horizon. 



Esh, GeV 
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,eV^ 









0.01 


0.005 


0.002 


I^M 










10 - 25 


1.47 


2.00 


2.02 


2.16 


25 - 63 


0.46 


0.88 


0.98 


0.64 


63 - 160 


0.14 


0.38 


0.24 


0.16 


> 160 


0.06 


0.09 


0.07 


0.06 


Vp. Vs 










10 - 25 


1.47 


1.14 


1.12 


1.08 


25 - 63 


0.46 


0.31 


0.30 


0.41 


63 - 160 


0.14 


0.09 


0.12 


0.14 


> 160 


0.06 


0.05 


0.06 


0.06 
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V. DISCUSSION AND CONCLUSIONS 

The rates of shower events are sufficiently high for a 
relatively crude detector as described in the Introduction. 
Four strings on radius of 15 m with length of 200 m each 
enclose a volume of 90 Kt of water or ice. This translates 
into a statistics of 190 events/yr above 10 GeV in the 
absence of oscillations and 270 events/yr for Am^ = 0.002 
eV^. Two years of measurement should be enough to 
observe the enhancement in the shower rate with good 
statistics if Am? is not lower than 2xl0~^ eV^. 

The requirements for the energy and angular resolu- 
tion of the detector are not very high. Table | gives the 
rates in wide bins to demonstrate the relatively low sensi- 
tivity to the energy resolution of the detector. Smearing 
of the rates shown in Fig. |^ with a Gaussian distribution 
mimicking energy resolution of 30% does not change sig- 
nificantly the detectability of Vr appearance. Restricting 
the measurement to showers closer to vertical direction 
enhances the appearance effect as the peaks and valleys 
in the Vt flux become more noticeable, although it de- 
creases the statistics of the signal events. An angular 
resolution of order 10° should be good enough for the 
measurement of the rates in tt steradian given in Table |. 

The presented calculation is not accurate enough to 
determine the actual detector rates because it assumes 
that: 

a) Vfj^{Vfj) CC interactions are identified and are not 
counted as shower events. In practice some of these 
events (at high y) would not be distinguishable and will 
increase the experimentally measured shower rate; 

b) all r decay energy not carried away by neutrinos con- 
tributes to Esh ■ Actually a fraction of the r decay energy 
is carried by muons and does not contribute to the shower 



Vu{^fi)vT-{vT) channel 



rate, especially in the /i^ 
with a branching ratio of 17.4% |11|. 
These effects are not included in the calculation because 
their strength depends very much on the detector design 
and can only be studied by detector Montecarlo codes 
that account for the Cherenkov light propagation in wa- 
ter or ice and detector efficiency. 

The absolute normalization of the predicted shower 
flux is not better than 20 - 25 % ||l|l. This does not, 
however, decrease significantly the sensitivity to oscilla- 
tion parameters because the energy spectra of the shower 
events are significantly different for Am? greater than 
about 2x10^'^ eV^. The spectral changes would be eas- 
ier to detect if the measurement could be extended to 
shower energies lower than 10 GeV, where the contribu- 
tion of oscillatated t neutrinos is not significant. 

The large uncertainty in the absolute normalization 
will not allow a definite conclusion for the case of oscil- 
lations in sterile neutrinos. A combination of the niuon 
neutrino disappearance with no increase of the shower 
rate will be, however, an important indication in favor of 
such oscillations. 

In conclusion, a measurement of the energy spectrum 



of neutrino induced showers will provide valuable com- 
plimentary information on the oscillations of atmospheric 
muon neutrinos. An increase of the shower rate would be 
a detection of the appearance of r neutrinos in v^^ — > v-r 
oscillations for Am^ greater than 2x10^'^ cV^. 

Acknowledgements The author acknowledges valuable 
discussions with R. Engel, T.K. Gaisser, F. Halzen and 
D. Seckel. This research is supported in part by U.S. 
Department of Energy contract DE FG02 OlER 4062. 



[2] 
[3] 
[4] 

[5] 
[6] 



[7] 
[8] 
[9] 



[10] 

[11] 
[12] 



Y. Fukuda et al, Phys. Rev. Lett. 81, 1562 (1998); 
Y. Fukuda et al, Phys. Lett. B433, 9 (1998); Y. Fukuda 
et al, Phys. Lett. B436, 33 (1998); Y. Fukuda et al, 
Phys. Rev. Lett. 82, 2644 (1999). 

K.S. Hirata et al. Phys. Lett. B280, 146 (1992); 
R. Becker-Szendy et al. Phys. Rev. D46, 3720 (1992). 
K. Scholberg, in Proc. of the 1999 Venice Workshop on 
Neutrino Telescopes (ed. M. Baldo-Ceolin) , in print. 
See a list of the proposed long baseline oscilla- 
tion experiments and a collection of documents at 



http:// www. hep. anl. gov/ndk/longbnews^ 



L. Pasquali & M.H. Reno, Phys. Rev. D59: 093003 
(1999). 

F. Halzen & J.G. Learned, Proc. Fifth International 
Symposium on Neutrino Telescopes, Venice (1993), ed. 
M. Baldo-Ceolin; see the contributions of C. Spiering, 
Z. Dzhilkibaev, L. Resvanis & L. Moscoso in Proc. Eight 
International Symposium on Neutrino Telescopes, Venice 
(1999), ed. M. Baldo-Ceolin, in print. 
L. Moscoso (for the Antares Collaboration), in Ref. Q. 
M. Cluck, E. Reya, and A. Vogt, Z.Phys. C67, 433 (1995) 
T. Sjostrand, Comp. Phys. Comm. 39, 347 (1986); ibid 
43, 367 (1987); T. Sjostrand, JETSET version 7.4, pro- 
gram manual. 

V. Agrawal, T.K. Caisser, P. Lipari & T. Stanev, Phys. 

Rev. D53, 1314 (1996). 

C. Caso et al. Eur. Phys. J. C3, 1 (1998). 

W. Frati, T.K. Gaisser, A.K. Mann & T. Stanev, Phys. 

Rev. D48, 1140 (1993) 



4 



